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The study of time-dependent, many-body transport phenomena is increasingly within reach of ultra-cold 
atom experiments. We show that the introduction of spatially inhomogeneous interactions, e.g., generated by 
optically-controlled collisions, induce a conducting-to-nonconducting transition in the transport of atoms in 
ID optical lattices. Specifically, we simulate the dynamics of interacting fermionic atoms via a micro-canonical 
transport formalism within both mean-field and other approximations. For weakly repulsive interactions, a quasi 
steady-state atomic current develops that is similar to the situation occurring for electronic systems subject to an 
^. external voltage bias. We also find that this atomic current is robust against the details of how the interaction is 

switched on. As the interactions exceed a threshold value, a mean-field conducting-to-nonconducting transition 
occurs due to energetic constraints. This transition is preceded by the atomic equivalent of negative differential 
conductivity observed in transport across solid-state structures, but here the eventual nonconducting state is one 
ON that emerges dynamically. 
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Figure 1. Schematic plot of the experimental set up we simulate: Non-interacting ultra-cold fermions are loaded into the lowest energy state 
of a ID homogeneous optical lattice. At t — 0, a focused laser beam then induces on-site repulsive interactions of the Hubbard type on the 
left half of the lattice. Due to the imbalance of interactions in the middle, a mass current is induced. We emphasize that there is no need to 
"mechanically" set the atoms out of equilibrium by tilting a potential nor is it necessary to introduce additional species of atoms or dissipation. 
Moreover, this setup may also be implemented to study quench dynamics 1131 . The gray dots emphasize that atoms may be in a superposition 
of different quantum states. 

I. INTRODUCTION 

Advances in experimental studies of quantum transport of ultra-cold atoms in optical lattices HHU draw attention to different 
aspects of non-equilibrium physics. In conventional condensed matter settings, an electrical current is induced by applying a 
voltage difference across the sample. The interactions of the charge carriers are usually homogeneous, except near impurities. 
In contrast, ultra-cold atoms are charge neutral and so far the motion of an atomic cloud has been generated by placing the cloud 
away from its equilibrium position [1] or by a sudden shift of the minimum of its trapping potential or its distortion [2— 4|. The 
interactions are generated by either adding another species of atoms 1 1 1 or tuning magnetic -field controlled collisions 13] 0). 

In this paper we instead suggest a method to drive a mass current of atoms via local tuning of interactions. This, in turn, reveals 
interesting phenomena otherwise difficult to observe in conventional solid state systems. A key to realize interaction-induced 
transport is controllable inhomogeneous interactions - a novel possibility offered by optically tunable collisions of ultra-cold 
atoms. The optical Feshbach resonance (OFR) [5-8 1 is a promising technique for controlling interactions locally using focused 
laser beams. For instance, the spatial modulation of density using the OFR of bosons [ 6 1 demonstrates the power of controllable 
inhomogeneous interactions. Other exotic equilibrium structures may also be generated by using this tool |9|. 

We also point out that these optically controlled interactions can be time dependent by turning on/off the induced interactions at 
certain times. The sudden imbalance of interactions within the cloud drives the system out of equilibrium. Thus, instead of using 
external driving mechanisms such as a voltage bias, it is possible, with suitable initial conditions and patterns of interactions, 
to induce a mass current. We show that for reasonably weak interactions the mass current of fermions is similar to a charge 
current in electronic systems. At a threshold interaction strength (dependent on the filling) a conducting-to-nonconducting 
transition occurs. This transition is due to a mismatch of energy spectra between the interacting and non-interacting parts of 
the system. Preceding the transition is a many-body counterpart of negative differential conductivity observed in solid-state 
structures ifTOlfTTI . where changes in the carrier density or subdivisions of Brillouin zones cause non-monotonic dependence of 
the current on the external field strength. Recently it has been shown that negative mobility of cold atoms could be demonstrated 
in temporally modulated optical lattices [12|. Here, instead we propose that one may observe negative differential conductance 
in a static optical lattice by tuning the interactions in real space. 



II. INTERACTION-INDUCED TRANSPORT 

In our approach - which explicitly follows the dynamics of mass transport - we implement the micro-canonical formalism 
(MCF) lfl4l [TBI as in Ref. Ifl6l [T71 which monitors the evolution of the correlation matrix. When applied to a system of 
noninteracting atoms on a ID optical lattice suddenly losing the atoms on the right half of the lattice, MCF shows very different 
dynamics for bosons and fermions Ifl6l . While the bosonic current decays to zero in the thermodynamic limit, the fermionic 
current exhibits a quasi steady-state current corresponding to a plateau in the current as a function of time. Importantly, the 
MCF is designed for finite closed quantum systems so it is particularly suitable for studying the dynamics of ultra-cold atoms. 
Although our proposed setup can be applied to bosons, to explore their dynamics theoretically one has to consider also the effect 
of the quasi-condensate which renders the analysis more involved. Therefore, in this paper we focus only on fermions. 

We consider N pa (er =\, I) two-component fermions, with N p ^ — N p ±, loaded into a ID optical lattice of size N, see Fig.fT] 
This type of finite-length ID lattice may be generated by inserting a thin optical barrier [ 18 1 into a ring of optical lattice [ 19 1. 
The resulting C-shaped lattice is geometrically identical to the setup here. This confining potential creates a uniform lattice and 
does not require a background harmonic trapping potential parallel to the lattice. The background interactions are assumed to 
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Figure 2. (Color online) (a) Nl,ct and (b) AN^ a — AA r | CT (t = 0) with TV = 128 and N pcr = 64 for each species from mean-field dynamics. 



be negligible. The system is initially described by the tight-binding tunneling Hamiltonian Hq = —tJ2 



(ij).a C ia^ja 



, where (ij), 



i, c i(T (Ci a ) denote nearest-neighbor pairs, the hopping coefficient, and the creation (annihilation) operator of site i, respectively. 
The unit of time is to = fi/t. 

The initial state corresponds to the lowest energy state of H , with no correlations between up and down spin states (see, 
e.g., Ref |[T6l \TT\ ). The system is then set out of equilibrium by introducing interactions among atoms on the left half of the 
lattice by, e.g., optically-induced collisions using a focused laser beam. Here we concentrate on moderate repulsive interactions 
and non-integer filling. This allows us to focus on the dynamics induced by the interactions and avoid unnecessary confusions 
about possible equilibrium phase transitions such as the BCS instability or the Mott-insulating phase in uniform and static 
interacting systems [20|. We model an instantaneous switch-on of the interactions with a sharp interface between the interacting 
and non-interacting regions and relax the former condition later on. 

The Hamiltonian generating the dynamics is 
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Here h in = c ia ^ a 



Ci a , U is the onsite repulsive coupling constant, L denotes the left half of the lattice, and a is the opposite of 
a. This model should be appropriate for moderate lattice depth ||2TI . The equations of motion for the correlation matrix are 
i(9(cj .c,- .) /dt) = ([cj CT , H e ]cj a ) + (c ia [cj a ,H e ]}, where [•, •] denotes the commutator of the corresponding operators. The 
explicit expression can be found using standard anti-commutation relations. After some algebra, one gets 
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(hi) and n ia = (h ia ). We solve Eq. |2]) both at the 
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mean-field level and by adding higher-order correlations. 

The mean-field level solution can be obtained by Wick decomposition of {c\ s CivC i(J Cj a ) as (c] 5 .Cjct-)(cJ .c,v) since no spin-flip 
mechanisms are present. This is the standard Hartree-Fock approximation. The equations are closed after this approximation. 
Thus the dynamics after the interaction is switched on can be monitored by integrating Eq. |2 1 with the initial condition c^ CT (i = 
0) given by the lowest energy state of the non-interacting Hamiltonian. We update Eq. ( 2| in a symmetric fashion so that 

( c L c j>> = (4a c js) at an y time- 
In addition to the current I a — 2tIm(c] v/2 .Cjv72+i,<t)> we aiso evaluate the particle number on the left half lattice, N^^ — 

(Nl<t) an d its fluctuations AJV£ CT = (iv£ ff ) — (TVt CT ) 2 , which is another property that can be determined experimentally. Here 

N La = YsicL ™io- Explicitly, 
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Right before the interactions are switched on, there can be initial number fluctuations due to the wave nature of the initial 
quantum state. We thus present AJV^,, — AN^,Jt = 0) which reflects the change of the number fluctuations due to the 
interaction-induced dynamics. 



III. RESULTS AND DISCUSSIONS 



Figure [2] shows the particle number and its fluctuations of the left half lattice with N = 128 and N pa — 64 for each species. 
As U/t increases, one can see that instead of inducing more fermions to move from the left to the right, above a certain threshold 
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Figure 3. (Color online) (a) I a and (b) (J CT ) of the fermionic case with TV = 128 and N pa = 64 for each species from mean-field dynamics. 
The inset of (b) shows G e ff = /~ 1 d(I< T }/dL r (in units of 2ir/h) for small U/t as a function of N pa for TV = 128, where / = TV pCT /TV is the 
initial filling factor. The dashed line shows Go ~ 0.152 in this unit. The energy spectra of the left (L) and right (R) half lattices when U is 
switched on are shown for U/t = 7 (c) and U/t = 8 (d). Here TV = 128 and N pa = 64. 



value the transport is no longer observed within a reasonable time scale (e.g., t < (N/2)to). This corresponds to a conducting- 
nonconducting transition and it can be found at other ratios of / = N p /N. In general the threshold U c /t decreases as / 
increases. We will show that this is because for higher filling the onsite interaction energy dominates more easily. One can see 
this transition even better in the average current itself. 

We note that it has been shown that with suitable initial conditions non-interacting fermions can develop a quasi steady-state 
current (QSSC) |fl*6ll22l . This is so also for the present interacting case. The fermionic current is shown in Figure |3| a) for 
TV = 128 and N pa = 64 for each species. The plateaus shown in I a {t) indicate the existence of QSSCs. To smooth over 
small oscillations, we define an averaged current as (I a ) = (l/30to) J 10t ° dtl a (t) and plot (I a ) as a function of U/t for 
N pa = 32, 64, 96 in Fig. plb). For small U/t we clearly see that the averaged current increases linearly with U/t. Then the 
dependence deviates from a linear form before the conducting-nonconducting transition. For even larger U/t, (I a ) decreases 
until no finite averaged current can be found. The non-monotonic dependence of (I a ) on U resembles the negative differential 
conductivity observed in conventional solid-state devices iTTOl ITTTl . 

The interaction term thus acts as a bias on the left half lattice if U/i is small. Here we investigate the validity of this 
statement. By extracting the slope of the linear part of (I a ) for small U/t, one can study how efficient the interaction can induce 
a current. The inset of Fig. Bib) shows the effective conductance G e // = f~ 1 d(I a )/dU as a function of TV pCT for a fixed TV, 
where / = N pa /N is the initial filling factor. As N ptT — > TV the system approaches a band insulator |23] and there is only a 
tiny range of U/t where a finite (I a ) exists. This could contribute to some inaccuracy of the extraction of the slope close to 
N p(7 — > TV. Nevertheless, as shown on the inset of Fig. 131b) the effective conductance is close to the ideal quantized conductance 
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Figure 4. (Color online) Currents induced by different ways of switching-on the interaction from mean-field dynamics: A sudden switch-on 
of the interaction (dashed line), linear switch-ons with switching times t m — 5io (black) and t m — 10io (red), and a multi-step switch-on 
(dot-dash line). The inset shows the time dependence of U for the four cases. Here TV = 128 and N p<7 — 64. 



Go = (2nh)~ l for arbitrary initial filling. This establishes the feasibility of using a weak inhomogeneous interaction as an 
efficient driving force for inducing a current and its correspondence with a bias, A/i = fU. 

In the weak interaction regime, we found another interesting phenomenon. By considering different ways of switching on the 
interaction, one may wonder whether the system will reach the same magnitude of the QSSC. We simulate different switching- 
on scenarios by introducing a time scale t m such that the interaction grows from zero to its full value during t m and remains 
a constant after that. Figure [4] shows that despite different time-dependences of the interaction, the system always reaches the 
same magnitude of the QSSC. We have also checked other more complicated functional forms for switching on the interaction 
and found the same QSSC. Importantly, even if the system is over-excited by a spike in the time dependence of U, when U 
comes back to a constant value the height of the QSSC coincides with the case of a sudden switch-on of the interaction as shown 
on Fig. [4] Thus the magnitude of the QSSC is robust (memory-less) against different ways the interaction is switched on. This 
feature makes the interaction-induced transport appealing for atomtronics since a steady output (the QSSC) is insensitive to the 
transient behavior of the driving force (the interaction). 

When the optically-induced interactions are switched on, it is possible that the onsite potential could be shifted due to the 
increase of the kinetic energy of the atoms interacting with the incoming photons. We study possible effects by adding an extra 
term V J2ieL a ^ic t0 H e of Eq. ([T| to simulate this effect with a positive V. This shift of the onsite potential acts like a bias 
so one might expect that more current will flow to the other side. However, the conducting-nonconducting transition is robust 
against this positive potential shift. The energetic mechanism described below will make it clear why a positive potential shift 
will not remove the transition, thus showing that it is not sensitive to this type of modification of the Hamiltonian. 

We now discuss the reason behind the transition, which is energetic in nature. In order to obtain a "macroscopic" steady-state 
current, aN particles must be transferred from the left to the right halves, where a is some small proportionality constant. 
However, this entails a change in energy, as the particle density decreases on the site with interactions, and this effect must be 
compensated by an energy change due to rearrangement of the single particle states. To be more concrete, immediately after 
interactions are turned on, the system is still in the ground state of H Q . However, the energy has shifted to 



E = E Q + U^(n ta n i9 ) = E Q + U^(n ta ) (n l9 ) , 



(4) 
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where Wick's theorem has been applied since the system is in the ground state of Hq. The subscript of (■ • • ) denotes the time 
when the expectation is taken. For half-filling, this gives an energy E = Eq + UN/8. For arbitrary filling /, it is 



E « E + fUN/2, 



(5) 



where the expression is approximate since not all site occupation numbers will be exactly / (for each spin component) . 

The energy at some later time will still be E (since no further changes have been made to the Hamiltonian), but will have the 
form 



E = (H ) t + UJ2(^)t(n l a)t 
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in the mean-field approximation. Assuming the aN atoms are taken uniformly from the left half, the energy is 

E « (H )t + f 2 UN/2 - afUN. (7) 

Taking the difference of Eq. (|7]i and Eq. (|5), and using energy conservation, gives (-ffo)t — Eq f=a afUN. However, the quantity 
(Ho) t — Eo can not be arbitrary, since we have a finite bandwidth. Regardless of whether one has an interacting or noninteracting 
state, the largest energy change due to the Ho contribution is 4taN, corresponding to taking the aN particles from the lowest 
part of the band (on the left half) to the highest part (on the right half). This gives 

afUN < AtaN — ► U < U/f (8) 

and thus results in a threshold value of U. Beyond this value of U, a macroscopic current can not flow. We remark that here we 
consider the lowest energy band of the optical lattice. If the interaction energy exceeds the band gap separating different bands, 
the system may re-enter a conducting state. 

To demonstrate this energetic mechanism, Fig. [3jc) and (d) show the energy spectra of the left (L) and right (R) half lattices 
when U is switched on for TV = 128 with N pa = 64. When U/t = 7 (panel (c)), there is still an overlap between the two spectra 
and exchanging atoms in an energy-conserving fashion is possible. For U/t = 8 (panel (d)), there is no overlap between the 
two spectra and the system evolves into a nonconducting state. This theory also explains why the nonconducting state is robust 
against an additional positive onsite bias. Any further increase of the onsite energy separates the two energy spectra of the left 
and the right sides even farther so the nonconducting state remains. Importantly, this energetic mechanism applies to other types 
of setups for studying transport phenomena. For example, by introducing a step-function bias to the lattice potential and tuning 
the bias, there is also a conducting-nonconducting transition of the same nature (see Appendix [A). 

We note, however, that other possibilities may occur in interacting systems for large U/t when a homogeneous conducting 
state is no longer favorable. From Eq. (l6j, states could be generated that have an oscillating density which would give enough 
"energetic relief" to allow a macroscopic current to flow [24|. In the mean-field solution, however, the state of the system 
is steered into a nonconducting state that is a dynamically generated zero current state. Since the conducting-nonconducting 
transition emerges dynamically, it should not be confused with conventional transitions of equilibrium states. For example, 
there are several differences when the dynamically-generated nonconducting state is compared to the Mott insulating state in 
equilibrium ll20l : i) the nonconducting state could occur at arbitrary filling while the Mott insulator needs integer filling; ii) the 
mechanism behind the nonconducting state is energy conservation while that behind the Mott insulator is the minimization of 
ground state energy; Hi) the nonconducting state depends on its initial state and evolution, while the Mott state is the ground 
state of the corresponding Hamiltonian. 

IV. HIGHER-ORDER CORRELATIONS 

We finally check whether the conducting-nonconducting transition survives when higher-order correlation functions (beyond 
using the Wick decomposition) are considered. The equations of motion of two-particle correlation functions were preserved and 
the three-particle correlation functions were decomposed into two- and single-particle functions as shown below. Importantly, 
we found that the conducting-nonconducting transition is still observable. We thus emphasize that even with the inclusion of 
higher-order corrections, the system dynamically evolves into the nonconducting state, thus supporting that this energetically 
driven transition should be observable experimentally [25|. 

To monitor the dynamics, we integrate the equations of motion while respecting the symmetry between the two spin indices. 
The equations of motion of the correlation matrix elements (or single-particle correlation functions) now include two-particle 
correlation functions. Above, we closed the set of equations using the Wick decomposition to write two-particle correlations 
as products of single-particle correlation functions, which is the Hartree-Fock (or mean-field) approximation. We now consider 
the case where the three-particle correlation functions are decomposed into products of single-particle correlation functions and 
two-particle correlation functions, e.g., 
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This allows us to include all of the two-particle correlation functions in the simulation. Since there are three ways to decompose 
the three-particle correlation functions, we chose to average the contributions of each in order to reduce the numerical instability 
of the simulation. 

We have computed the currents by solving the mean-field equations of motion and the equations that include two-particle 
correlations for N — 40 and N pa = 20. The results are shown in FigureBland they agree qualitatively. More importantly, while 
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Figure 5. (Color online) The currents from the mean-field approximation (dashed lines) and from the simulations including two-particle 
correlations (solid lines). Black (red) lines correspond to U/t — 2 (20). Here N = 40 and N pa = 20. 



finite smooth currents are found in both approaches for small U/t, the current oscillates rapidly with a small amplitude around 
zero for large U/t. Figure p] shows the average current including higher-order calculations, and is similar to the mean field case 
in Figure p[b): There is an increase then a decrease in (I a ) as U/i is increased, which is also observed for other higher order 
approximations. There are differences, however, with the mean field results, as shown in FigurelT] especially within the center of 
the lattice. Nevertheless, there is no indication that a conducting state emerges for large U (e.g., a state with density oscillations 
sufficient to counter balance the energetic constraints discussed above). 

We note that inclusion of the two-particle correlations is only a short time approximation and the simulations can become 
unstable for intermediate to long times. We have thus restricted these simulations to short times and are therefore limited to 
averaging over shorter time scales than above. This means that the averaging partially truncates oscillations and the effect of 
this partial truncation is to give a finite value to (I a ) even when a longer average would give zero. This will be more prominent 
for shorter averages, i.e., in the absence of transient effects, one expects that the boundary will give a residual contribution of 
(I a ) oc l/U, as U gives the frequency scale of oscillations (for large U). The current versus time shows that this is indeed what 
is happening: There are faster oscillations as U is increased (and also as / is increased). When averaging to get (I a ), one will 
not obtain zero for large U but rather some residual effect because of the short time average. The inset of Figure [6] shows this 
tail. For the filling greater than 1/2, the average current is consistent with this interpretation, and thus the nonconducting state 
is forming. For filling equal to 1/2, the results are ambiguous, as the oscillations are both slower than the higher filling cases 
and many-body correlations are more important. Nevertheless, for / = 1/2, the average current still decays to zero, which is 
indicative of the nonconducting state. 

Furthermore, our results suggest that for small U/i the transient time to reach a quasi steady state is only a few t while the 
duration of the QSSC can extend much further than that. Thus for reasonably large lattices, the observation of an interaction- 
induced QSSC and the conducting-nonconducting transition is feasible. 



V. CONCLUSION 



In summary, we have shown that one may use time-dependent inhomogeneous interactions of ultra-cold atoms to explore non- 
equilibrium physics not easily realizable in conventional solid-state setups. In particular, we find a conducting-nonconducting 
transition, where the nonconducting state emerges dynamically. Importantly, the current in the weakly-interacting regime is 
robust against possible transient behavior of the driving force. Our micro-canonical simulations may be compared with other 
numerical methods such as the time-dependent density-matrix renormalization-group method [26 1 to further investigate the 
correlation effects. Our work not only sheds light on fundamental physics out of equilibrium, it may also lead to different 
designs of devices in the thriving field of atomtronics J27ll . 

CCC acknowledges the support of the U. S. DOE through the LANL/LDRD Program. MD acknowledges support from the 
DOE grant DE-FG02-05ER46204 and UC Laboratories. 




Figure 6. (Color online) (I„) of the fermionic case with N — 64 and N pa = 32, 48, 56, 60 from the simulations with two-particle correlations. 
Similar to Figure [3] the average current initially increases with U and then decreases. Here, the averaged current is computed from i = to 
t — 2.5to. The inset of the plot shows the average current plotted on a log scale, which shows that it does indeed vary as 1/U a where a tends 
towards 1 as the filling increases (a w 2.95, 1.59, 1.43, 1.29 for N pa = 32, 48, 56, 60 respectively). Partially truncated oscillations are thus 
likely responsible for the residual value of (I a ) at large U. 
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Figure 7. (Color online) (a) Difference between the onsite correlations (hi a ) as calculated using the mean field approximation and using the 
two-particle correlations and (b) the relative difference between the onsite density-density correlations {hi a fiig). While the two approximations 
can differ in the value of these correlations, there is no evidence that density oscillations (beyond those seen in the mean field) emerge from 
the simulations. Thus, the formation of the nonconducting state for large U is expected. 



Appendix A: Comparison of bias vs. interaction 



When a bias is applied to a single-species noninteracting fermions in a uniform ID lattice, the Hamiltonian is Hb 
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bias (fis = 0) and then it evolves according to Hb- For TV = 128, we define (I a ) = (l/30trj) J wt ° dtl a (t), where to = H/t. 
Figurep|(a) shows (/„•) as a function of fis- One can see that above /zb/£ = 4 no finite averaged current can be found. There is 
a conducting-nonconducting transition and the upper limit is independent of N p /N, which can be seen by a similar analysis to 
Eqs. (4)-(8) in the main text. 

As with the interaction case, energy conservation prohibits particles on the left with energy higher than the highest energy 
level of the right from entering the right half of the lattice. This occurs when fis is larger than the bandwidth of the original 
uniform lattice, which is At. The particle distributions of the two half lattices when the interaction is switched on can be found as 
follows. The initial correlation matrix Cij(t = 0) of the whole lattice is X)p=i Uo(hP)Uo(p,j), where Uq is the unitary matrix 
that diagonalizes the hopping Hamiltonian. After finding Ul and Ur that diagonalize the Hamiltonian of the left and the right 
half lattices with eigenvalues El and Er, the particle distributions are given by n L / R (k) = J2i j^l l «(^ l I R)kjUik<kj(t = 0), 




Figure 8. (a) (I a ) as a function of /j,b for N = 128 and N p = 32, 64, 96. The curves for N p = 32, 96 coincide due to the particle-hole 
symmetry. Inset: The currents for /is/t = 1, 2, 4 (labeled next to each curve) with N p — 64. (b) Particle distributions of the left half (black) 
and the right half (red) for /is/t = 3 (solid line) and 4 (dashed lines). The symbols show energy levels. Here N = 128 and N p — 64. (c) 
Same as (b) but with N p = 96. 



where k = 1, • • ■ , N/2 labels the energy levels E L / R {k). Fig. [81(b) and (c) show the particle distributions for the left and right 
half lattices for /.is/i = 3 and 4 with N p = 64 and 96. One can see that when the bottom of the left half energy band moves 
above the top of the band of the right, no current is expected since energy is conserved. 

Taking the same definition for (I a ) for the interaction-induced case, we plot (I a ) as a function of U in Fig. [9] There is also a 
conducting-nonconducting transition, but the threshold value depends on the filling N pa /N . The left half lattice has Hamiltonian 
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To find the energy levels of Hl at t = 0, we replace hi„ by 



its expectation value, which is known from Cij t0 -(t = 0) and we assume the particle number of the two species are equal. After 
diagonalizing H^ = Y1uj)£L a c la- c jc + U^2 ieL ni S hi a , we found Ej J (k) with k = 1, • • • , k and the unitary matrix Ul- 
The right half lattice remains noninteracting and one can find En(k) and Ur. The distribution functions are then given by 
n L /n{k) = T,i,jeL/R( U L/R)kjUikCij(t = 0). 

Fig. [91(b) and (c) show the distributions of the right and left half lattices for N pa = 64 and 96 and N = 128. The nonconduct- 
ing states emerge slightly above U/t = 7 (U/t = 5) for N pa = 64 (N pa — 96). The distributions for U/t below the threshold 
value still show observable overlaps between the filled states of the left half lattice and the empty states of the right half lattice. 
When U/t is above the threshold value, the distributions of the two half lattices are separated so energy conservation forbids 
transport. We thus observe nonconducting states in both bias-induced and interaction-induced transport due to energy-level 
mismatches of the two half lattices. In the interacting case, however, there may exist states that could carry a current, but our 
simulations show that the nonconducting state emerges dynamically in the setup we examine. 
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Figure 9. (a) (J) as a function of (7 for JV = 128 and N p — 32, 64, 96 from the mean-field approximation, (b) Particle distributions of the 
left half (black) and the right half (red) for U/t = 7 (solid line) and 8 (dashed line). The symbols show energy levels. Here TV = 128 and 
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